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Abstract: Cooperative communication in wireless sensor network (WSN) explores the 
energy efficient wireless communication schemes between multiple sensors and data 
gathering node (DGN) by exploiting multiple input multiple output (MIMO) and multiple 
input single output (MISO) configurations. In this paper, an energy efficient cooperative 
MIMO (C-MIMO) technique is proposed where low density parity check (LDPC) code is 
used as an error correcting code. The rate of LDPC code is varied by varying the length 
of message and parity bits. Simulation results show that the cooperative communication 
scheme outperforms SISO scheme in the presence of LDPC code. LDPC codes with different 
code rates are compared using bit error rate (BER) analysis. BER is also analyzed under 
different Nakagami fading scenario. Energy efficiencies are compared for different targeted 
probability of bit error p b . It is observed that C-MIMO performs more efficiently when 
the targeted p b is smaller. Also the lower encoding rate for LDPC code offers better error 
characteristics. 

Keywords: cooperative technique; LDPC; BER; MIMO; wireless sensor networks 
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1. Introduction 

Recent advances in micro-electro-mechanical systems technology have enabled the development of 
wireless sensor nodes in a wireless sensor network (WSN). These tiny sensor nodes are able to sense, 
process and communicate with each other [1,2]. Since the battery capacity in each node is limited and 
the goal is to maximize the lifetime of the network, there are strict energy consumption constraints in 
WSNs [3]. The size of sensors is typically small but the functions inside the sensor are complex. Recent 
hardware advancements allow more signal processing functionality to be integrated into a single sensor 
chip. RF transceiver, A/D and D/A converters, base band processors, and other application interfaces 
are integrated into a single device to be used as a smart wireless node. A wireless sensor network 
typically consists of a large number of sensor nodes distributed over a certain region. Monitoring 
node (MN) monitors its surrounding area, gathers application-specific information, and transmits the 
collected data to a data gathering node (DGN) or a gateway. Energy issues are more critical in the 
case of MNs rather than in the case of DGNs since MNs are remotely deployed and it is not easy to 
frequently change the energy sources. Therefore, the MNs have been the principal design issue for energy 
limited wireless sensor network design. One prospective solution is the use of MEMO [4,5] for energy 
efficient design with a targeted probability of bit error at the receiver. Also LDPC-coded MIMO optical 
communication is mentioned in [6]. But the MIMO techniques require complex transceiver circuitry 
and signal processing leading to large power consumptions at the circuit level. Moreover, physical 
implementation of multiple antennas at a small-size sensor node may not be feasible. The solution came 
in the form of cooperative MIMO (C-MIMO) [4-8]. C-MIMO is a kind of MIMO technique where 
the multiple inputs and outputs are formed via cooperation in a network of single antenna nodes. The 
sensors cooperate with each other to form a MIMO structure and in fact lead to better energy efficiency 
and smaller end-to-end delay. The basic idea of C-MIMO was first proposed by S. Cui in [4]. Later this 
idea has been improved in [5] by Jayaweera considering channel estimation (training overhead) in the 
DGN side and is further modified in [9] by Y. Gai and in [8] by M. Rakibul. 

The issue of applying error control codes in WSNs is the topic of recent interest. The performance of 
block codes and Viterbi decoded convolutional codes is investigated in [10,11]. The iterative decoding 
algorithm using turbo code is used to prolong the network lifetime [12]. Low-density parity-check 
(LDPC) codes are more reliable than the block and convolutional codes and are serious competitors of 
turbo codes. In particular, LDPC codes exhibit an asymptotically better performance than turbo codes 
and admit a wide range of trade-offs between performance and decoding complexity [13]. Sartipi and 
Fekri [14] compare the performance of the LDPC codes and the Reed Solomon (RS) codes [15]. From 
the recent works, it is known that LDPC codes are attractive in WSNs because of their applications in 
compression, joint sourcechannel coding and distributed source coding [14,16]. However, to the best 
knowledge of the authors, there has been no document on the implementation of LDPC encoder/decoder 
in a wireless sensor node using cooperative communication. More precisely, none of the recent works 
have addressed the problem of reducing the energy consumption using error control coding. In this 
paper, LDPC code is incorporated in cooperative communication as an error control code. Later the idea 
is compared with SISO communication. In Section 2, the system model is shown and the error correction 
using LDPC code is analyzed in Section 3. Section 4 shows the energy model for both cooperative MIMO 
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and SISO considering error correction codes. In Section 5, simulation results are shown and discussed. 
Finally, Section 6 concludes the paper. 

2. Cooperative MIMO Communication 

2.1. System Model 

The system model for C-MIMO communication is a centralized wireless sensor network where there 
is a data gathering node (DGN) and several clusters with several sensors in each cluster. Sensors in 
one cluster transmit the data to the sensors in adjacent cluster and step by step the data reach the DGN. 
Figure 1 shows the cluster to cluster communication where two clusters are shown. The system considers 
N t number of sensors in the transmitting cluster, N r number of sensors in the receiving cluster and one 
antenna is placed at one sensor. Also, each element in the channel matrix H is assumed to be a zero-mean 
circularly symmetric complex Gaussian random variable with unit variance and can be considered as 
follows. 

hi Nr ^ 

h2N r 

\ h^ti h^ t 2 . . . hx t N r J 
Figure 1. System model for cluster to cluster communication in wireless sensor network. 




Cluster 1 Cluster 2 



The problem here is stated from the receiver point of view, so a loss model is used to estimate the 
received energy. To calculate the total energy consumption, both the circuit and transmitter power are 
taken into count. The same transmitter and receiver blocks shown in [4] are used in this paper. Source 
coding, pulse shaping, and modulation block are as well omitted from the design. Throughout the 
paper, a system with narrowband, frequency-flat Rayleigh fading channels and perfectly synchronized 
transmission/reception between wireless sensor nodes is assumed. The fading is assumed constant during 
the transmission of each frame. In our model, a sensor with high residual energy is deployed as a cluster 
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head and it remains the cluster head until the network dies. The cluster head broadcasts its status to 
the other sensors in the network. Each sensor node determines to which cluster it wants to belong by 
choosing the cluster head that requires the minimum communication energy. Once all the nodes are 
organized into clusters, each cluster head creates a schedule for the nodes in its cluster. This allows the 
radio components of each non-cluster-head node to be turned off at all times except during its transmit 
time, thus minimizing the energy dissipated in the individual sensors. 

2.2. Cooperative Communication 

The physical phenomena monitored by sensor networks, e.g., forest temperature, water 
contamination, usually yield sensed data that are strongly correlated. Data aggregation is the tool by 
which the correlated data size can be significantly reduced depending on the correlation factor. Figure 2 
explains the cooperative communication where the sensors at cluster 1 send the information data to the 
cluster head os cluster 2. At the first step, the sensors at cluster 1 send the data to their cluster head. 
The cluster head then aggregates the data in the second step. After the aggregation, the cluster head 
send the aggregated data back to all the sensors in that cluster. This is the step three in cooperative 
communication. At this stage, all the sensors at cluster 1 have the same information data. At the fourth 
step, the sensors transmit the aggregated data to the cluster 2. After receiving the data at the receiving 
cluster, sensors at cluster 2 transmit the received data to their cluster head locally and complete the 
cooperative communication. 

Figure 2. Cooperative communication. 
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3. Error Correction Codes in Wireless Sensor Network 

Error control coding (ECC) introduces redundancy into an information sequence u of length k by the 
addition of extra parity bits. Several different types of ECC exist, but we may loosely categorize them 
into two divisions: (1) block codes, which are of a fixed length n c , with n c — k parity bits, and are 
decoded into one block or codeword at a time; (2) convolutional codes, which, for a rate k/nc code, 
input k bits and output n c bits at each time interval, but are decoded in a continuous stream of length 
L » nc- Block codes include repetition codes, Hamming codes [17], Reed Solomon codes, and BCH 
codes [18]. Short block codes like Hamming codes can be decoded by syndrome decoding or maximum 
likelihood (ML) decoding by either decoding to the nearest codeword or decoding on a trellis with the 
Viterbi algorithm or maximum a posteriori (MAP) decoding with the BCJR algorithm. Algebraic codes 
such as Reed Solomon and BCH codes are decoded with a complex polynomial solver to determine 
the error locations. Convolutional codes are decoded on a trellis using either Viterbi decoding, MAP 
decoding, or sequential decoding. 

Another categorization is based on the decoding algorithms: (1) non-iterative decoding algorithms, 
such as syndrome decoding for block codes or maximum likelihood (ML) nearest codeword decoding 
for short block codes, algebraic decoding for Reed Solomon and BCH codes, and Viterbi decoding or 
sequential decoding for convolutional codes; (2) iterative decoding algorithms, such as turbo decoding 
with component MAP decoders for each component code, and the sum product algorithm (SPA) [19] 
or its lower complexity approximation, min-sum decoding [20], for low density parity check codes 
(LDPCs). The non-iterative decoding category may be further divided into hard and soft decision 
decoders; hard decision decoders output a final decision on the most likely codeword, while soft 
decision decoders provide soft information in the form of probabilities or log-likelihood ratios (LLRs) 
on the individual codeword bits. Viterbi decoding can be either hard decision or soft decision, with 
a 2dB gain in performance for soft decision decoding. Category (2) are all soft decision algorithms 
by nature, as iterative decoding requires soft information as a priori input for each iteration. Iterative 
decoding algorithms provide significant coding gain, at the cost of greater decoding complexity and 
power consumption. With the recent technological advancements, all these ECC techniques can be used 
in WSN. However, LDPC code is considered in this paper as an ECC tool at WSN for its superior error 
correcting capabilities. 

3.1. Low Density Parity Check Codes 

Low density parity check codes are codes specified by a matrix containing mostly 0's and relatively 
few l's. A standard bipartite graph based ensemble which is shown in [21,22] is used in this paper. The 
code length is designated by n and number of constraints by m. Therefore, there are n variable nodes 
and m check nodes. Each variable node corresponds to one bit of the codeword and each check node 
corresponds to one parity check equation. Edges in the graph connect variable nodes to check nodes and 
represents the nonzero entries in H matrix. The term "low density" conveys the fact that the fraction of 
nonzero entries in H is small, in particular it is linear in the block length n, as compared to "random" 
linear codes for which the expected fraction of ones grows like n 2 [23]. 
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For regular codes, the corresponding H matrix has 8 r ones in each row and 8 C ones in each column. 
It means that every codeword bit participates in exactly 5 C parity-check equations and that every such 
check equation involves exactly S r codeword bits. Low density parity check codes have been constructed 
mostly using regular random bipartite graphs. 

Example 1. Here is an example of a regular parity check matrix with S r = 6 and 5 C = 3. 

"111001100010" 

111110000001 
H _ 000001110111 

100100011101 

010110111000 
_001011001110_ 

The bipartite graph corresponding to this parity check matrix is shown in Figure 3. 

Figure 3. Bipartite graph corresponding to a regular parity check matrix. 



Variable nodes 




Check nodes 



For irregular codes, 6 r and 5 C are not fixed for every row and column of the parity check matrix. We 
consider that the irregular bipartite graph has a maximum variable side degree 5 r and a maximum check 
side degree S c . 

Example 2. The following H matrix is an example of an irregular parity check matrix with a 
maximum 5 r = 6 and a maximum 5 C = 3. 

"101001100010" 

110110000001 
H _ 000001110111 

100100011101 

010110111000 
_001011000100_ 

The bipartite graph corresponding to this parity check matrix is shown in Figure 4. 

For this paper, H matrix shown in Figure 5 is used for simulation. This H matrix is a special matrix 
used in 802.1 In standard. WSN is energy constraint in nature and the sensors work as intermediate 
devices when the data are transferred from a designated area to the data gathering node (DGN). Since 
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decoding can be performed in the DGN, energy efficient decoding technique is not a concern for this 
paper. Encoding is one critical issue considered in the wireless sensor network. In this paper, Richardson 
encoding scheme is used as a tool for using LDPC code in WSN and is explained in the next subsection. 

Figure 4. Bipartite graph corresponding to an irregular parity check matrix. 

Variable nodes 



Check nodes 
Figure 5. H matrix used in the simulation. 
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3.2. Richardson Scheme as the Encoding Technique 

The encoding method proposed by Richardson et al. [13] assumes that H can be converted to an 
approximate lower triangular matrix. The authors worked on an m x n parity check matrix H over F 
where n is the number of variable nodes and m is the number of check nodes. The parity check matrix 
H is transformed in the form of 



\ 



\ 



H 



A B T 
C D E 



(1) 



where A is (m — g) x (n — m), B is (m — g) X g, T is (m — g) X (m — g), C is g x (n — m), D is 
g x g, and E is g x (m — g) where g is denoted as the gap. Further, all these matrices are sparse and T 
is lower triangular with ones along the diagonal. This matrix is multiplied from the left by 
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/ 0 

-ET' 1 I 



(2) 



And the H matrix is found as 



A B T 

-ET~ X A + C -ET~ X B + D 0 



(3) 



They then break the codeword as x = (s, pi,P2) where s denotes the systematic part, p± and p 2 denote 
the parity part, p\ has length g, and p 2 has length (m — g). After that, the equation Hx T = 0 T is used to 
state the following two equations 



As T + Bp\ + Tj% = 0, 



(4) 



(-ET^A + C)s T + (-ET^B + = 0. 
Taking {—ET^B + D) as nonsingular, it is concluded that 

pi = —(—ET~ l B + D)- 1 {-ET~ 1 A + C)s T . 



(5) 



(6) 



P-2 



-T-\As T + BpT 



(7) 



By using the step by step procedure, it is shown that the complexity of calculating p 1 and p 2 are 
0(n + g 2 ) and 0(n) respectively. The matrix used in our simulation can also be written in ALT form 
and is shown in Figure 5. 

4. Energy Model for Cooperative Communication Using LDPC Code 

The total power consumption P T for a single node consists of two main parts, namely, the power 
consumption of all the power amplifiers Ppa which is a function of transmission power P out , and the 
power consumption of all other circuit blocks Pc- Thus one can write 

P T = P PA + P c . (8) 
The power consumption of all the power amplifiers can be calculated using the following equation 

P PA = (l + a)P out , (9) 

where a — — 1), where r/ is the drain efficiency [24] and £ is the peak to average ratio. When the 
channel only experiences a k th power path loss with additive white Gaussian noise (AWGN), P out can 
be calculated using the link budget relationship as follows. 



out 



E b R b x ±-L—MtN,, 



(10) 
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where E h is the average energy per bit required for a given bit error rate (BER) specification, R b is the 
transmission bit rate, d is the transmission distance, G t and G r are the transmitter and receiver antenna 
gains respectively, A is the carrier wavelength, Mi is the link margin compensating the hardware process 
variations and other background noise, Nf is the receiver noise figure defined as Nf = where N r is the 
power spectral density (PSD) of the total effective noise at the receiver input and No is the single-sided 
thermal noise PSD at the room temperature. 

The power consumption in the circuit block includes transmitter and receiver power consumption 
P ct and P cr , respectively. This power consumption is due to several power blocks such as P m ix, P syn , 
Pfut, Pfur, Plna, Pi fa, Pdac, an d Padc which are the power consumption values of the mixer, the 
frequency synthesizer, the active filters at the transmitter and at the receiver side, the low noise amplifier, 
the intermediate frequency amplifier, the D/A and A/D converter, respectively. The power consumption 
block for error correction is not considered as it is same for cooperative case and SISO case. The total 
energy consumption per bit can be written as 

„ _ (Ppa + Pc) nu 
k, ht 5 , (11) 

K b 

where R b is the actual bit rate and can be replaced by R b * = f ~p Nt R b when pN T training symbols 
are inserted in each block to estimate the channel at the receiving cluster or DGN side. The block size 
is equal to F symbols and can be obtained by setting F = \T c Rs\ , where Rs is the symbol rate and 
Tn is the fading coherence time. The fading coherence time can be estimated from Tn = , 1 3 7= where 
the maximum doppler shift f m is given by f m = v with v being the velocity and A being the carrier 
wavelength [25]. The total energy consumption is estimated by multiplying E bt by the number of bits L 
to be transmitted. Now we develop the mathematical model where we estimate total energy consumption 
for cooperative communication. 

The total energy consumption in cooperative case is modeled as 

N t -l N t N t j 

EcO = ^iPi + Pda + ^enc ^ — 7i 

i=l i=X i=l 

N t 



r 
i=i 



Nt L 

i=l 

N t N r 



i=l i=l 



(12) 



The energy per bit E\ is needed to transmit the data from sensors to the cluster head. Ed a is the 
energy dissipation per bit required in the cluster head for data aggregation. It depends on the algorithm 
complexity and can be expressed as 

Co + Ci x L for O n 
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where L is the number of transmission bits and C 0 , C\ and C 2 are coefficients depending on the software 
and CPU parameters. E enc is the encoding energy per bit and is taken 1 fxJ [26]. Ef denotes the local 
transmission energy cost per bit for transferring the aggregated data to the remaining active sensors, 7 is 
the percentage of remaining data after aggregation and it reflects the correlation between data amongst 
different sensors, r is the rate of LDPC encoding. Since the use of a rate r = 1/2 makes the size of 
the data after encoding, 2 times the original data size, the ^ term is used to represent the data size after 
encoding a message size of Lj with rate r. The same energy per bit Ef is needed to transmit a command 
signal from the cluster head to the selected sensors. After receiving all the bits, the nodes encode the 
transmission sequence according to some diversity scheme, such as the STBC. E l M denotes the energy 
cost per bit for the long haul MIMO transmission [4]. J2f=i "7"7« i s divided by the optimal bit size of the 
long haul transmission b m i mo to find the number of symbols present in the received signal. The number 
of symbols is then multiplied by the optimal bit size of the local transmission h ir to find the total bit 
length. Ej is the energy per bit required to transmit the data from a sensor to the cluster head at the 
receiver side. iV r is the number of sensors at the receiving cluster. 

For the SISO approach, sensors transmit their data to the cluster head and as there is no burden for 
channel estimation, the cluster head will transmit all the aggregated data directly to the destination node 
without any cooperation. So the total energy consumption becomes 

N t -1 N t N t j N t j 

Esiso = LiE i + Eda Y. u+ E ™ c S 7^ +E sJ2 -7^' (14) 

i=l i=l i=l i=l 

where E l s denotes the SISO long haul transmission and can be calculated as a special case of MIMO 
transmission with Nt — 1 and iV r = 1. In both SISO and C-MIMO case, fixed constellation size is used. 
Since the encoding energy using Richardson scheme is same for both C-MIMO and SISO approach, it 
is not considered in the equation for C-MIMO and SISO. 

5. Simulation Results and Discussion 

In order to get the total communication energy consumption, the average energy per bit required for a 
given BER P b , E b needs to be determined. In this approach, the value of E b is found by using a numerical 
search. Ten thousand randomly generated channel samples are taken and averaged to find the desired bit 
error rate at each transmission distance. The value of the constellation size is kept fixed. For the long 
haul communication, SISO is taken as a special case of MIMO structure. A list of system parameters 
used in our simulation is shown in Table 1 where the power consumption values of various circuit blocks 
are quoted from [4] . 

5.1. Energy Issue 

Total energy consumption and energy efficiency are the key terms to evaluate the energy efficient 
performance. For simulation, it is considered that all the sensors in a cluster are transmitting the same 
data size of Li = 10 kb. The simulation is performed based on the cluster size of Nt = 4. In Figure 6, the 
total energy consumption over distance is shown for cluster to cluster data transmission. From Figure 6 
it is clear that the cooperative MIMO is more energy efficient than SISO transmission. The simulation 
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is taken for two different code rates r = 1/2 and r = 3/4. When the code rate increases, parity bits 
compared to message bits are reduced. Therefore, total energy consumption reduces. This is verified in 
Figure 6. 

Table 1. System Parameters. 



f c = 2.5 GHz 
G t G r = 5 dBi 
B = 10 KHz 

p b = io- 3 

N f = 10 dB 
P syn = 50.0 mW 
E da = 5 nJ/bit/signals 
P filt = 2.5 mW 



rj = 0.35 

iV 0 = -171 dBm/Hz 
k = 2 for local com 
k = 3 for long haul com 
Mi = 40 dB 
P mix = 30.3 mW 
P LNA = 20 mW 
P fUr = 2.5 mW 



Figure 6. Total energy consumption over distance. 
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5.2. Delay Issue 

The total delay required is defined as the total transmission delay. For a fixed transmission bandwidth 
B, we assume that the symbol period is approximately T$ ~ 1/B. The total delays in the case of SISO 
communication is defined as 

( Nt L 1 Nt L \ 
Tsiso = T s V ^ + V - 7i + *da, (15) 

Where 6^ is the transmission bit size at the transmitter side local communication and bsiso is me 
transmission bit size for long haul SISO transmission. td a is the time taken for data aggregation. 
The total delays in the case of cooperative MIMO communication is defined as 



Sensors 2011, 11 



9898 



N t 



Nt Li 



-CO 



=1 Lt i=\ lt 




>MIMO 



+7 V1Et 



(16) 



Where and t^ a are the channel estimation and data aggregation delays respectively. The term 
T s J2f=i TT~ x% f° r me delay due to the local transmission from sensors to the cluster head. The next term 
is due to the local transmission from cluster head to the sensors. Tqt — - — Y\ *, — term is caused by 

a b MIMO L^ii=\ r J 

the long haul MIMO transmission. The next term is due to the local transmission at the receiver side. 
The assisting nodes first quantize each symbol they receive into n r bits, then transform all the bits into 
symbols using bi r and transmit to the cluster head to do the joint detection. 

The delay difference is calculated using the following equation. We assume the value of t ch w 0. 

DD = Tsiso — Tco 

N t r N t Li 



bsiso ~ r 
1 Nt T 

t — y. —it 

Omimo , r 

1=1 



E 



i=l 



Tc 



i Nt 



(17) 



The value of n r is chosen at the receiver based on the optimized transmitted constellation size. The 
delay difference is a measure of delay performance by which the cooperative MIMO can be compared 
with SISO. Positive delay difference indicates the SISO is facing larger delay compared to C-MIMO. In 
Figure 7, delay difference is compared where proposed C-MIMO outperforms SISO after 60 meters. 



Figure 7. Delay difference over distance for code rate r 
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5.3. Constellation Size Issue 



Since Energy efficiency = {E S iso — Ec-mimo} I 'Esiso, positive energy efficiency indicates 
{E S iso > Ec-mimo}- In Figure 8, energy efficiency is simulated over distance for different 
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constellation sizes. The simulation results show that for rate r = 1/2, cooperative MIMO outperforms 
SISO after 80 meters for constellation size b = 1 whereas it takes 10 meters for constellation size b = 8. 

Figure 8. Energy efficiency for different constellation sizes. 




5.4. Bit Error Rate Issue 



Using the parity check matrix provided in IEEE 802. lira standard shown in Figure 5, comparative 
error performance studies have been taken for different encoding rates and are shown in Figure 9. Also 
the C-MIMO is compared with SISO in the \ rate case. The codeword length is kept fixed and the 
number of decoder iteration is taken as 100. Bit error rate (BER) is taken as performance parameter in 
this paper. BPSK modulation and AWGN channel are used for the simulation. Like the other wireless 
channels, simulation using cooperative MIMO shows similar outcomes. 

Figure 9. Bit error rate over SNR curve for different encoding rate. 
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The same BER analysis is taken in Figure 10 in a Nakagami fading channel scenario. The result 
shows that the decrease in Nakagami coefficient m degrades the error performance. 

Figure 10. Error performance for different Nakagami coefficient m in Nakagami fading 
channel. 
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Figure 11. Energy efficiency for different targeted probability of bit error p b . 
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5.5. Reception Quality Issue 

Targeted BER is the parameter that indicates the reception quality of the signal. The cooperative 
MEMO communication used in this paper is simulated with a fixed value of targeted BER. Figure 1 1 
shows that the change in targeted BER changes the efficiency of cooperative communication over SISO 
transmission. Result shows that the cooperative communication is more energy efficient than SISO 
transmission in smaller targeted BER. 

6. Conclusions 

Energy efficient data transmission is one of the key factors for energy constraint wireless sensor 
network. An energy efficient cooperative technique considering low density parity check codes is 
modeled and simulated using Matlab. The results show that the cooperative communication outperforms 
SISO transmission at the presence of error correction code. The energy efficiency remains almost 
unchanged in different encoding rates but it largely varies with the change in constellation size. BER 
analysis is also taken to show the similar error characteristics in the cooperative MIMO environment. 
Data with smaller encoding rate shows better BER results than larger encoding rate for a fixed SNR. 
Simulation is also performed in the situation of a fading environment. It is also found that cooperative 
communication is more energy efficient than SISO transmission in smaller targeted BER. Therefore it 
can be concluded that cooperative MIMO with LDPC can be a good choice for high reception quality 
signals. 
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